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Abstract 

In many extensions of the Standard Model, finite temperature computations are com- 
plicated by a hierarchy of zero temperature mass scales, in addition to the usual thermal 
mass scales. We extend the standard thermal resummations to such a situation, and 
discuss the 2-loop computations of the Higgs effective potential, and an effective 3d 
field theory for the electroweak phase transition, without carrying out high or low 
temperature expansions for the heavy masses. We also estimate the accuracy of the 
temperature expansions previously used for the MSSM electroweak phase transition in 
the presence of a heavy left-handed stop. We find that the low temperature limit of 
dealing with the left-handed stop is accurate up to surprisingly high temperatures. 
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1 Introduction 



Electroweak baryogenesis in the Minimal Supersymmetric Standard Model (MSSM) is 
a viable option for explaining the matter-antimatter asymmetry observed in the present 
Universe, provided that there is a mild hierarchy between the right and left-handed 
stop masses 0-0. The dominantly right-handed stop should be lighter than the top 
in order to make a strong transition, yet the left-handed stop should be rather heavy, 
~ 1 TeV, in order to raise the Higgs mass upper bound to ~ 110 GeV. Various details 
of the electroweak phase transition in this regime are constantly being investigated [0] . 

Here we will be concerned with the thermodynamics of the phase transition. In the 
perturbative approach, this problem is approached by computing the effective potential 
for the Higgs field to some order in the loop expansion. In general, such a computation 
in a weakly coupled gauge theory faces two problems: (i) The system has a hierarchy of 
mass scales (27rT, gT, g 2 T), which spoils a straightforward perturbative computation. 
Historically, this was observed by finding large "linear terms" at 2- loop level ||, which 
were then shown to be absent after an appropriate resummation ||. (ii) At momenta 
of the order of the lowest of the mass scales (g 2 T), the system is also inherently non- 



perturbative [10 



The resummations needed at 2-loop level for dealing with the heavy scale 2itT were 
discussed in detail by Arnold and Espinosa fll |. However, the problem can also be 



dealt with in another way, namely by constructing a sequence of effective field theories 
by integrating out, to a given order in perturbation theory, the scales 2nT, gT |TJ], [L3|]. 
This construction is highly accurate in the Standard Model [|TJ], TZ\. The final theory 
is three-dimensional (3d), purely bosonic, and contains only the momentum scale g 2 T. 
A perturbative analysis of the 3d theory automatically reproduces the results of the 
resummed 4d effective potential, but the theory can also be studied efficiently with 



relatively simple lattice simulations IfLql , to account for the non-perturbative part. 

The problem we consider here is the observation that the hierarchy of mass scales can 
be even more severe in extensions of the Standard Model such as the MSSM. Indeed, 
there one tends to have new mass parameters that are not related to the temperature 
in the same way as ma is in the Standard Model, where m# ~ gT c . In particular, 
as mentioned above, one prefers rather large left-handed squark mass parameters, say 
niq ~ 1 TeV. Previously, the effects of tuq have been considered (on the 2-loop or non- 
perturbative level) only in the high temperature expansion, or in the extreme limit 
rriQ ^> 2nT where the finite temperature effects decouple completely. 

Our objective here is to treat in some detail the general situation tuq ~ 2nT. First 
of all, we discuss how the resummations used previously need to be changed in such 
a situation (Sec. |3|). We then show with a simple example how the full resummed 
2-loop effective potential could be computed without any temperature expansions re- 
lated to tbq, and how the result can be used for a 2-loop computation of the mass 
parameter of an effective 3d field theory (Sec. f|). Finally we consider a particular 
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observable sensitive to rriQ, the critical temperature of the electroweak phase transi- 
tion, and estimate the accuracy of the high and low temperature expansions employed 
earlier on (Sec. [5]). We conclude in Sec. ^|and discuss several possible extensions of the 
computations presented in this paper. The expressions used for the 1-loop tadpole and 
bubble, as well as 2-loop sunset graphs are discussed in the appendices. 



2 Parametric conventions 

In order to be explicit yet concise, we illustrate the situation with a simple model 
reminiscent of the scalar sector of the MSSM. We take 

C = m] i H^H + m 2 u U* a U a + m 2 Q QlQ a 

+ hjH^HU* a U a + h 2 2 H^Q a QiH + h 3 (AH^Q a U a + R.c. )+.... (2.1) 

Here H is an SU(2) doublet, U an SU(3) (anti-)triplet, while Q changes under both 
groups. We ignore gauge interactions for the moment. We assume that h\ ~ /12 ~ 
h 3 ~ g are small couplings, and m^,m^ ~ (gT) 2 . It is also important to specify the 



order of magnitude of the dimensionful parameter A in Eq. ( |2.1| ). In this paper we 
work under the assumption that 



\A 



2 



lAf^^-g 2 , (2.2) 



2 

m 



which simplifies the procedure considerably. 

In the imaginary time formalism, the fields in Eq. ( |2.1| ) can be divided into Matsubara 
modes. We assume that the only light modes are the zero Matsubara modes of H, U . 
The non-zero Matsubara modes of H, U have effectively a mass parameter > (2irT) 2 . 
For the field Q, we assume that ttiq itself is large, mg ~ 2tcT, so that even the zero 
Matsubara mode is heavy. If mg ~ gT, then the zero Matsubara mode of Q is light 
as well and the procedure is the one described in [jXTJ. If ttiq ~ 2nT/g, on the other 
hand, Q can be integrated out at T = with exponentially small corrections. 

The issue of resummation can now be formulated as follows. Due to the presence of 
the heavy mass scales, the n = modes of H, U can receive radiative corrections as large 
as the tree-level terms, 5m 2 HU ~ g 2 Tnq + g 2 T 2 . Such corrections have to be resummed. 
In fact, close to the phase transition point, the effective mass parameters m 2 HeS Ueff can 
be even smaller, of the non-perturbative magnitude ~ (g 2 T) 2 . Then resummation has 
to be extended to the 2-loop level. Non-zero Matsubara modes, or the field Q, on the 
other hand, do not require resummation ||11|| , since the mass corrections g 2 T 2 ,g 



' 2 m 2 n 



are according to our convention small compared with the tree-level terms. 

This statement can be formulated more precisely as follows. Let us write down the 
effective Lagrangian obtained after integrating out all the heavy modes. The light fields 
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Figure 1: The 2-loop graphs considered. The blob means a counterterm. Different 
particles are denoted by Pj = H,U,Q. 

being the n = modes of H, U, the form of the Lagrangian is 

£ eff = m 2 HcS H^H + m 2 UeS U* a U a + h 2 lcS H^HU* a U a + ... . (2.3) 
Our aim is now to compute expressions of the forrrf] 

m HeS = m H + Cig 2 m 2 H + c 2 g 2 m 2 u + c 3 g 2 m 2 Q (l + c 4 g 2 ) + c 5 g 2 T 2 (l + c 6 g 2 ), (2.4) 
where the Cj-'s are numerical coefficients. 

3 Leading order resummation 

In order to carry out the resummation explicitly, let us consider the effective potential 
of the theory. To illustrate the procedure, it is enough to consider only the field H, 
keeping the expectation value of U at zero. Introducing (H) = (0 , 0) T / v2 changes 
the mass spectrum of the system, m = m + d^m. We are interested in a certain 

range of (f) ~ 0...T. Then, in the case of heavy modes, m, and we can expand 

in 5^771, while in the case of light modes we cannot. For the purpose of illustration, let 
us suppress m for the light modes here. 

Then, in the standard thermal case, the 1-loop and 2-loop contributions to the 
effective potential behave at small <ft a s 

~ T 2 (6m^ 2 + T(5m (j> ) 3 + ... J (3.1) 
SV 2 . loop ~ g 2 THm <t) + g 2 T 2 (5m 4> ) 2 + .... (3.2) 

The statement of resummation is now that the dominant 2-loop terms, the "linear" ones 
~ g 2 T 3 Sm ( p, arise from a badly convergent series which can be resummed into a better 
convergent one 0. The way the resummation proceeds is obvious from Eqs. (|3.1|) , 



(|3.2| ): the non-analytic 1-loop and 2-loop terms combine to 

T{5m^f + g 2 T 3 6 m4> -> T{g 2 T 2 + 5m 2 ) 3/2 . (3.3) 



3 In a gauge theory there are also corrections of order g 3 T 2 
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This corresponds simply to the corrections of order g 2 T 2 in Eq. (|2.4j) . The extension 
we make here is that when mg ~ 2nT, the contribution to be resummed goes from 
g 2 T 2 to a non-trivial function g 2 T 2 f(rriQ/T, \A\ 2 ). 

In order to proceed systematically, we write the mass parameters related to the light 
modes as 

™H = ™HeS ~ °rm H , TTljj = m UeS - « r %, (3.4) 

where ^# e fr; m f/ e ff appear in the propagators, and 5 r m 2 H) b r rn\j are treated as interac- 
tions. Denoting the heavy modes by solid lines and the light modes by dashed lines, 
the graphs 

..Q. --0-- 

suggest that 

S r m 2 H = 3hll nj ta(mu) + 3h 2 2 I{m Q ) + 3hl\A\ 2 I{m Q ) - J„^ (0) , (3.5) 
b T m 2 v = 2h 2 1 I n ^ {m H ) + 2hl\A\ 2 \l{m Q )-I n ^{0)}. (3.6) 



Here I,I n ^o are tadpole integrals defined in Eqs. ( |A.1| ), ( |A.11| ), and we have made 



use of m^,m^ <C m 2 Q. Note that the fact that \A\ 2 is small, Eq. (|2.2| ), implies that 
wave function corrections need not be considered, since their effect would be of order 
~ h\\A\ 2 m 2 H ~ g 4, m 2 H , beyond Eq. (|2.4j) . For the same reason, we have dropped any 
w? Hl m^j dependence in the terms proportional to |A| 2 . 

In addition to the mass parameters of the scalar fields, resummation of course also 
affects the zero components of the gauge fields. In fact, as is well known || ||, in 
the Standard Model the latter effect is more important for physical observables such 
as the strength of the phase transition, while the former is important particularly for 
the critical temperature. We do not discuss infrared dominated observables such as 
the strength of the phase transition, nor gauge fields, to any length in this paper, but 
let us nevertheless note that the contributions of H, U, Q to the Debye masses of the 
SU(2) and SU(3) fields A , C are, in the presence of trq ~ 2irT, 

$rm 2 Ao = g 2 T-^(l n ^ {m H ) +3/ T (m Q )), (3.7) 

5 rm 2 Co = ^T^(^oK) + 2/ T (m Q )), (3.8) 

where gs is the SU(3) gauge coupling and It is defined in Eq. (|A.6| ). In addition to these 
terms, the Debye masses of course contain the usual gauge and fermion contributions. 

In order to now show that the procedure introduced in Eqs. ( |3.4| ), (|3.5|), ( |3.6| ) is a 
consistent one, we need to demonstrate that all "linear terms" at 2-loop level cancel, 
and the remainder is quadratic in ^m^. Recalling that we have set the quartic Higgs 
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self-coupling to zero (at tree-level) for the purpose of simplicity, we get for the shifts 
in the mass parameters (Qi(2) denote the upper (lower) SU(2) component of Q) 



S<t>m 2 HeS 



0. 
0. 



(3.9) 
(3.10) 



Note that due to the assumption |v4| 2 ~ g 2 we can ignore all corrections involving A 
here, since the corresponding 2-loop contributions are at most of order ~ h 2 h\\A\ 2 ~ g % . 
We will denote (m^ efr ) 2 = m 2 HcS + S4>m 2 HeS , etc. 

Linear terms in the effective potential arise from graphs of the types (H), (U), (HU), 
(HQ), (HUQ) in the notation of Fig. [3]. Denoting by 7 3d the 3d tadpole in Eq. fl7TT2; ) 
and by H the bosonic sunset integral in Eq. (|C.1|) , we obtain 



(H) + (U) = -25 r m 2 ¥ J 3d (m^ eff ) - 35 r m 2 / / 3d (mJ eff ), 
(HU) + (HQ) = Qhlnmi^HmU) + 3h 2 2 I(mU) [/(mjj + I(m% 2 



(HUQ) 



-3hi\A\ 



H(™.Q1, m HeS, m Uefi) + #0™Q2> m HeK, m Ues) 



(3.11) 
(3.12) 

(3.13) 



We then expand these contributions in S^m. Employing the expansions 
I(m 



1 



m? 



(-/(m Q ) + /^ (0) 



I(m Q ) - S^m 2 Q D(m Q ) + 0(S^m Q )\ (3.14) 
J ( m fs) = 4^o(ra e ff) + ^(rafff) - ^rn,e ff £) n ^o(^eff) + 0(W 4 , (3.15) 

J 3d("4cff) + J 3d(«4eff)) 

+/3d(m^ efr )/ 3 d(^ eff )] + 0(^m) 2 , (3.16) 
where D, D nj L are from Eqs. ( |B.5| ), ( |B.10| ) and we have used Eq. (|C.17| ), we find that: 

• there are linear terms oc hd(^Hes) ^3d( m Ues) which however all get cancelled, 
when the choice in Eqs. (|3.5|), ( p.6[ ) is made for 5 r m 2 H , b r m\;. 

• there is an infrared sensitive contribution, quadratic in the masses, from the 
Matsubara zero modes in the graphs (HU), (HUQ): 



(HU) + (HUQ)\ 



IR 



6(hl - h 2 3 \A\ 2 )I 3d (mU)h d (mi eS ). (3.17) 



The appearance of h\ — h\\A\ 2 corresponds to coupling constant resummation 
which we however do not discuss in any detail here, since the corresponding 
effects are in principle beyond the accuracy of Eq. ( ^.4| ). Similarly, the graph 
(HUQ) also produces terms of order ~ /i 4 |A| 2 2 , again beyond Eq. 
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• finally, there are ultraviolet sensitive (not from the zero modes) quadratic terms 
from the graphs (HU), (HQ): 

(HU) + (HQ)\ VY = -6hlS^mu e gInM m H) D nM m u) 

-3h 2 2 S^m 2 Q I n ^(m H )D(m Q ) + 0(5^m)\ (3.18) 

To summarize, we have observed that all linear terms get cancelled when the thermal 
counterterms are chosen according to Eq. ( |3.11| ). The remainder involves quadratic 



terms, which can either come from the ultraviolet or the infrared. 



4 Next-to-leading order 

We next evaluate the 2-loop contributions from the remaining graphs, and expand them 
again in 5^m\ however, these graphs do not involve contributions linear in 8$m. The 
graphs left are the sunsets (HQQ), (HUU), as well as the 1-loop graphs (H), (U), (Q), 
where the blobs are now the bilinears obtained from the coupling constant counterterms 
after the shift of H []. 

After an expansion in 8^m, we obtain 



(H) + (U) + (Q) 

(HQQ) 
(HUU) 



1 1 



2 (4tt) s 



9(/4 + h A 2 )I n ^(m H ) + Qh\l nj LQ(mu) 



+12/4/(m Q )] + 0(5^mf, 
-3h\^H(m Q , m Q , 0) + 0(8^mf 

~h\4> 2 H(m* H%s , mfj eS , mfj eS ). 



(4.1) 
(4.2) 

(4.3) 



The numerical expression of H(rriQ, tuq, 0) is discussed in appendix C. As to the graph 
(HUU), on the other hand, we recall that it arises completely from the zero Matsubara 
modes (H = H 3d + 0(m/T)), and is thus a purely IR quantity |2T| . 

Adding all terms together from Eqs. (|3.18|) , ((O]), fl4~2^ ) and using Eqs. (|A.11|) , ( p3|) , 
( [B.lOp , ( |C.20| ), we obtain the 2-loop ultraviolet contribution to the 3d mass parameter, 



xuv 2 _-4 

2-loop " L Hefi— lh 2 



-6H vstc (m Q ,m Q ,0) + 12 I vac (m Q ) 



+ 



T 2 



(4tt) 2 
3 



3 1 
11 



(4*0 

- In ^5 — 3( In — + c 

4 /ij, V /i 



-hi 



In 



( m Q\ 



2 r , 

Q 1 



In 



+ 2 



■ (4.4) 



4 Mass counterterms do not contribute at the present order; terms proportional to tuq in them 
would, had we included self-interactions of the type ~ (WH) 2 , (U*U a ) 2 in Eq. (2.1). 
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Here the first line is a 2-loop vacuum renormalization correction of order g mn 

1 



ln T + W) ~ 2ie] k _0 - 34872274 ' (4 - 5) 



fl T = 4vre- 7B T « 7.0555T, l , } 

and Ti, V, 7i are functions defined in Eqs. ( |A.7| ), ( |B.8j ), ( |U.2 

On the other hand, the IR sensitive part of the effective potential is, from Eqs. (|3.17|) , 

®> 

SocpV = 6(/i? - /#4Vsd(^ (4.6) 

The divergence in if 3d (Eq. ( jC.2Q ) cancels against that from Eq. (|4.4|), m 2 HeS (f) 2 /2. 

Including also the 1-loop terms in Eq. (|3.5| ), we can now write down the complete 
mass parameter m 2 HeS with accuracy g^rrig, g A T 2 . In order to do so, let us first note 
that 1-loop radiative corrections generate couplings other than those in Eq. (|2.1| ), viz. 



5C = hlH^HQiQ a + X(H^H) 2 + . .., (4.7) 

which we have to include in the discussion for a moment. The corresponding contribu- 
tion in Eq. ( |3.5| ) is 5 r m 2 H = 6h 2 I(mQ) + Q\In^o(jn 2 H) ■ Furthermore, in order to cancel 
spurious /i-dependences, we should express the MS parameters in terms of physical ob- 
servables as in ||12j| . In this paper we will not consider actual physical pole masses etc, 
but simply some finite physical scale independent parameters () p h y s which, dropping 
all terms beyond the accuracy of Eq. fl2.4p, we define through the following relations: 



3 



2 



m H (ji) = m Hphys + j—^ [Kmjj + (h 2phys + h 3 \A\ )m Q \ [\n — T + 1 

+h 2 6H vac (m Q ,m Q ,0) - 12D vac (m Q )I vac (m Q ) (4.8) 



him) = ^W + ^iT^p 111 ^' ^(M) = h 2 2phys + hl-^\n-2, (4.9) 
h 2 M = h'+hi-^ln-^r, A(/2) = \ phys + hh 4 1 + hj) — ln -4,(4.10) 



where H v£LC (m Q , m Q , 0), D v£LC (m Q ), I vac {m Q ), are defined in Eqs. ( |U.19| ), ([BT6|) , ( |7H>1 ). 
Moreover, let us now declare h\ phys) A p h ys ~ 0. We then obtain the final expression for 
the 2-loop effective (bare) mass parameter m 2 HeS in the theory of Eq. (|2.1D: 



2 2 ^ 7 2 2 ( \ 



+ T 2 {\(h 2 lphys - h\\A\ 2 ) + 3 -(hj phys + hWM 2 )!^) 
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Figure 2: Left: The coefficient of T 2 from Eq. ( [4.11] ), along with the high-T and low-T 
limits. The parameters have been chosen as explained in Sec. |5]. Right: The ratio of 
the high-T and low-T limits to the non-expanded result. 



(4tt) 



rrin 4 e 



4 



3T 

3(hi— + c 
v V 



-^(l2X 1 (^) + ip(^) + 6W(^) 



(4.11] 



5 High-T and low-T expansions 



We now wish to employ Eq. ( |4. 1 1| ) to estimate in a non-trivial physical context the 
accuracy of the high and low temperature expansions in tuq/T. We can do this by 
inspecting the critical temperature T c of the phase transition. Let us recall that the 
leading (and next-to-leading in a gauge theory) terms in T c are perturbative 



thus ultraviolet dominated and particularly sensitive to tuq/T. Most of the physical 
characteristics of the phase transition, on the contrary, are infrared dominated and less 
sensitive to ttiq/T. We may also remind that in the MSSM the determination of T c is 
physically more important than in the Standard Model, since one has to address the 
question of whether other phase transitions could take place before the electroweak 
one, in particular a transition to the dangerous ?7-direction 

The transition will take place when m 2 HeS (fi = g 2 T) = #g A T 2 , where $ is some non- 
perturbative coefficient, to be determined with lattice simulations. We shall keep the 
physical parameters () p h y s fixed and vary tuq/T. It is then clear that the perturbative 
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contribution to T c can equivalently be inspected by considering the finite part of the 
coefficient of T 2 in Eq. (|4.11|) . We choose for simplicity /ii p h ys , ^2 P h ys ~ 1, h\ \ A | 2 ~ 0. 
The high and low temperature limits of 2i,V,H are given in Eqs. ( |A.9| ), (|B.9|) , ( |C.22|) . 
To be in accordance with the limiting procedures usually applied in the literature, we 
keep in the high temperature expansions terms up to logarithmic order, whereas in the 
low temperature expansions we simply replace the exponentially small corrections in 
Eqs. flOD , (gj , ( ET2p with zero. 

The numerically evaluated full expression for the coefficient of T 2 , as well as a com- 
parison of the high and low temperature expanded versions thereof with the full result, 
are plotted in Fig. |2|. We observe that the high temperature expansion gives typi- 
cally too large a coefficient of the T 2 -term, leading to too small a T c []. With the 
low-temperature expansion, on the other hand, T c is slightly too large. Furthermore, 
we observe that while naively one might have expected the crossover between the high 
and low temperature regimes to be close to the first non-zero Matsubara frequency 
at irtQ ~ 2ttT, the low temperature expansion is in fact perfectly sufficient already at 
rriQ > 3T, which is the case for realistic values tuq > 300 GeV. The fact that the high 
temperature expansion converges relatively poorly as early as at ttiq/T ~ 2 is due 
particularly to the 2-loop function 7i, whose behaviour is shown in Fig. El 



6 Conclusions 

In this paper, we have pointed out that standard thermal resummations should be 
extended in two different ways, when one goes from the Standard Model to a general 
MSSM. First of all, the left-handed stop irtq is typically of the order of magnitude 
~ 2nT. Then it cannot, a priori, be treated either in the high or in the low temperature 
expansion, but a more general function appears. Second, the presence of dimensionful 
trilinear couplings leads to the emergence of new "linear" terms coming from the scalar 
sunset diagrams. The results for the scalar thermal counterterms including both of 



these effects in the model of Eq. ( |2.1|) are shown in Eqs. (|3.5| ), (|3.6|) , while the scalar 
contributions to the Debye masses are shown in Eqs. ( |3.7| ), ( |3.8| ). In an effective theory 
approach such as the one followed in [|l9j (in it was assumed that rriQ ^> 2ttT, but 
the procedure can be extended to rriQ ~ 2ttT in a straightforward way), all these effects 
of course arise automatically, whereas in a direct computation of the 2-loop effective 
potential they should be explicitly taken into account. 

In the framework of effective field theories, we have also extended the resummation 
for the Higgs mass parameter to the next order beyond the effects described above. The 
mass parameter thus determined, including corrections of order ~ g 4 T 2 , could be used 



Numerically the relative effect is larger here than in the realistic MSSM, since in that case there are 
other terms in the coefficient of T 2 (such as gauge bosons) for which the high temperature expansion 
should work perfectly. 
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for a precise estimation of the critical temperature of the corresponding electroweak 
phase transition using 3d lattice Monte Carlo simulations. Let us stress that the 
only change with respect to previous effective 3d theories is in the expressions for the 
effective parameters, not in the functional form of the theory. 

Using these results, we have estimated the accuracy of the high and low temperature 
expansions used previously in the literature. Inspection of the critical temperature 
suggests that the low temperature expansion, whereby all finite temperature contribu- 
tions from heavy particles are simply left out, works well already at m > 3T for bosonic 
particles. Thus it should be completely clear that for the values tuq ~ 1 TeV of interest 
for obtaining a strong phase transition with experimentally allowed Higgs masses in 
the MSSM, the Q-field can simply be left out in all finite temperature contributions. 

The present results could clearly be extended in many directions. First of all, the 
restricted model we have employed here can be extended to the full MSSM with gauge 
fields and fermions in a straightforward way. Second, we have shown that the evaluation 
of the integrals appearing in the perturbative 2-loop effective potential is numerically 
feasible without any further temperature expansions — thus the complete 2-loop po- 
tential of the MSSM could in principle be computed, extending thus the results of f|, 
@~0i 0- Third, we have here considered explicitly only the effects of a heavy rug, 
while in the MSSM many other mass parameters could be heavy as well. In particular, 
Mi-, A* related to the gaugino and Higgsino mass matrices and also relevant for provid- 
ing sources of CP violation can have values for which neither high nor low temperature 
expansions are applicable. In the effective theory approach M2, fi can be easily in- 
cluded at 1-loop level without any temperature expansions [fL9| |, but this could now be 
extended to the 2-loop level. The accuracy of previous approximations with respect 
to the contributions from the second Higgs doublet with txia ft 100 GeV could also be 
explicitly checked. Finally, we assumed that the trilinear couplings are not exceedingly 
large, \A\ 2 <g 2 , an assumption which could be relaxed. 

We believe that as long as the existence of a Higgs particle with MSSM type cou- 
plings lighter than about 110 GeV is not experimentally excluded, these are worthwhile 
questions to consider. 
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A The tadpole 



For completeness, let us review here some properties of the tadpole integral, 

d 3 ~ 2t p 1 



Urn) = T 



(2nf^p 2 + p 2 + m 2 ' 



(A.l) 



where po = p b = 2imT for bosons, po — Pf = 7cT(2n +1) for fermions. We will need 
two types of subdivisions of I(m). In the first case, relevant for all fermions and heavy 
bosons, we write I(m) = J vac (m) + /^(m), where It vanishes at T = 0. In the second 
case, relevant for light bosons (to 2 ~ (gT) 2 ), we separate the contribution from the 
Matsubara zero mode into /3d (to), writing I(m) = J 3 d (to) + I n ^{m). We denote 



u p ,= (p 2 + m 2 y/ 2 , u p =(p 2 + (m/T) 2 y/ 2 , 
d 3 ~ 2e p n b{f) (uj Pii ) 



lT,b(f)( m i) 



(2k) 



3-2e 



U). 



p,i 



When the superscript () b j is left out from J, we assume the bosonic case. 

A heavy mass in the loop. Writing I b (m) = I vac (m) + Jt(to), we get 

,2 n 



/vac (TO) 

I T (m) 



= -ix 



-2c 



TO 



1=4 + 1)- 



-/i" 2e T 2 



(An) 2 \ e to 2 

l + e ( 2 -21n2 + ln^)]l 1 (^)- e l 2 ( 7 , ;i , 



-4 / <w 



,n 6 (u) f 



1 f 00 

— / dp p 2 In p 

ir 2 Jo 



.n b (u p ) 



The limiting values are 



Jl(?/) = 6-2^ + 8^( 1 + 21n 7~ 27 ^' 

zb(v) 



2/»i /ye 



,-2/ 



27T 7T ' 



3 L 7£ + C(2) 



(A.2) 
(A.3) 
(A.4) 



(A.5) 
(A.6) 
(A.7) 
(A.8) 



(A.9) 



-In i/, 

7T 



»>i [ye 



2n 7r 



(lny + 2- 7f ;-ln2), (A.10) 



where lE = 0.57721566, C'(2)/C(2) = -0.56996099. 
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A light mass in the loop. Writing I(m) = I nj Ao(m) + /3d (m), we obtain 



InM m ) 
hd(m) 



/i_2£ 12 (1 + ae 



T 2 



m 



^ (4vr) 2 \r 



- + ln^) + 0(e 2 ,em 2 ,m 4 ),(A.ll) 



-2c 



mT 



Air 



Ti 

1 + e(lnAr + 2 -21n2 



0(6 2 ). 



Here \1t is from Eq. (|4.5|) and [TT 



l^ + 27^-21n2-2^| 



(A.12) 



(A.13) 



The fermionic tadpole. Using the standard trick, the fermionic tadpole can be 
expressed in terms of the bosonic tadpole: 



I f (m) = I vac (m) + I TJ (m), (A.14) 
I TJ (m) = 2/ w (m) -J Tib (m) = 2- 1+2e J Ti6 (2m) -J T , 6 (m). (A.15) 

Defining T 1 j,X 2 j as in Eqs. ( |A.7|) , (|A.8|) but with rif instead of n b , we obtain 



r.'i 



fva 



tr.j ( m ) ^/r 2e T 2 ^ - [l + e (2 - 2 In 2 + In ^)] J lf/ (^) + eZ 2J (- 



(A.16) 



The high and low temperature expansions of X^j,X 2 j can be obtained from those of 
Xi,X 2 by noting that 



^ , ... , 1 /2m\ , „_ /2m\ 
X 2 (-j--X 2 (-)+ln2Z 1 (— ). 



(A.17) 
(A.18) 



B The bubble 



Let us then consider the 1-loop "bubble" diagram with two propagators, 

1 1 1 



%)K,m 2 ) 



p b(f) P 2 + m 2 P 2 + m 2 m 2 — m 2 



h(f){m 2 ) - 4m (mi) , (B.l 



where Pb(f) = (Pb(f),p) and we have taken the external momentum to zero. We can 
again write 



%)K,%) = D vac (m 1 ,m 2 ) + D Tmb ( f) (mi,m 2 ), 
L) vac (mi,m 2 



I + In ^ + ! _ m l + m l ln nil 



%(/)K,m 2 ) 



(47r) 2 \e m\m 2 

_ 2e T 2 1 
& ^ 2 m 2 - m 2 

12 



m 2 — m 2 m 2 



(B.2) 
(B.3) 

0(e). (B.4) 



The special case mi = m-i gives the derivative of I(m) with respect to m 2 : 

dl(m) 



D(m) 
P> wac {m) 

D T {m) 

\t 



dm 2 



D vac (m) + D T (m), 



(4tt) 2 Ve 



m z 



-2c 



(4tt) 2 



4 r dv m 

Jo u n 



with 



„<i 2vr j/ 
^ (y) = — + 21n— + 2 7i; , 
y 47T 



S/>1 



y 



We also need the derivative of I n ^o(m) with respect to m 2 : 



C The sunset 



dm 2 



K 2 



(4tt) 2 V6 ' ^/Z 2 . 



+ ln^ +0(e,m 2 ). 



Let us then consider the bosonic and fermionic 2-loop sunset diagrams 

1 1 1 



H b ( f) (m 1 ,m 2 ,m 3 ) = E E 

JPu f \JQ 



% U) T QbU) P 2 + m 2 Q 2 + m 2 (P + Q) 2 + m\ ' 



In the limit rrii/T -C 1, it is known that [11, 21 

T 2 (\ 



H b (m 1 ,m 2 ,m 3 ) 



-Ac 



+ In 



(47r) 2 \4e mi + m 2 + m 3 2 
Hf(m 1 ,m 2 ,m 3 ) = 0(m,iT). 



(B.5) 
(B.6) 

(B.7) 

(B.8) 

(B.9) 
(B.10) 



(C.l) 



+ -)+0(miT), (C.2) 
(C.3) 



Our objective here is to compute these diagrams in the case of general m^. We are 



aware of previous results in this direction in [22, 23 



General case. The method we employ for evaluating Hf,, Hf follows the standard 



procedure (see, e.g., [|22], |23|, |24|). The twofold sum over the Matsubara modes is first 
written as a threefold sum with a Kronecker delta, and the delta is then written as 
$(Po) = T Jo dxexp(ip x). The sums can now be performed, 



7 ^ , 

Pb(f) Pb(f) + u i 



_ n h{f) {ui) 



2uJi 



e ( ± } e 



(C.4) 
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The integral over x is then very simple. The outcome can be organized in a transparent 
form, when different types of contributions are identified with known expressions; the 
same result could also have been obtained from the rules of the real time formalism, as 



noted for the 3- loop bosonic basketball diagram in ||25|| . In the remaining integral over 
the spatial vectors p, q, we can perform at least the integration over z = p ■ q/(|p| |q|), 
leaving for numerics at most a rapidly convergent 2d integral over p 
Let us denote 



Pi Q 



q 



/ M K,m 2 ;m 3 ) 



rf 4-2ep 



21 ' 



In 



(27r) 4 - 2e [P 2 + m\\ [(P + Q) 2 + m. 
4(p 2 + m\)(q 2 + ml) — [m\ + m\ — m\ — 2pq) A 



4(p 2 + m 2 )(q 2 + m 2 ) — (m 2 + m\ — in 2 + 2pq) 2 
The explicit expression for n(Q 2 ; 777 2 , m 2 ,), often denoted by Bq, is well known: 



(C.5) 
(C.6) 



H(Q 2 ; m\, ml 
F E (Q 2 ;ml,m\ 



// 



-2e rl 



(4tt) : 



+ In 

e m\m 2 



+ 1- 



mf — 777,2 m 2 



, (C.7) 



777? + 777n , 777i 777 

1 + * I In — + — 



777^ 



777T 



m\ 



777 2 



Q 2 



111 



777! 
777 2 



+7)2 V ( mi + rn 2 ) 2 +Q 2 \l{m 1 - m 2 ) 2 + Q 2 hi 



(mi-m2) 2 +Q 2 
(mi+m 2 ) 2 +Q 2 



(mi— m 2 ) 2 +Q 2 
(mi+m 2 ) 2 +Q 2 



(C.J 



The absolute value inside the logarithm in j VA in Eq. ( |C.6| ) means that we take the 
real part of the expression; the imaginary part would anyway cancel against Im II. 
With this notation, we obtain 

H b (m lt m 2 , 777 3 ) = F vac (mi, m 2 , tt7 3 ) 

+ E ^T,f,(^)Ren(-777 2 ;777 2 ,777^) 



+ E 



-4c 



327T 4 Jo 



dpp 



aqq ■ ■ j P , q {mi, mj, m^, 



UJ, 



(C.9) 



H f (m 1 ,m 2 ,m 3 ) = H vac (m 1 ,m 2 ,m 3 ) 

+I T , b (m 3 )Re Ilf-m 2 .; mj, mj) - Irjim^Re n(-m 2 ; m 2 , m\) 



-k 





~32tt 4 

E 



3 /"DO 7 

cipp / dgg- 



-k 



32vr 4 Jo 



dpp / dqq 



f P ,q(m 1 ,m 2 ;m 3 ) 
fp^ q (m 3 ,m i ;m j ), 



(CIO) 
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where Ei^k = E(ij,fc)=(i,2,3),(2,3 1 i),(3 1 i,2) ) = E(ij)=(i,2),(2,i), and the zero tempera- 
ture contribution is 



-"vac 

mi, m 2 , m 3 J 



1 



(27r) 4 " 2e J (2vr) 4 - 2e P 2 + m\Q 2 + m\ [P + Q) 2 + m| 



(C.ll) 



The only ultraviolet divergences are in if vac (mi, m 2 , m 3 ), which has 1/e 2 , 1/e poles, and 
in the IPs, which have the pole /i~ 2e /(167r 2 e), cf. Eq. ( |C.7| ). 

Let us also mention a few words about the numerical evaluation of the 2d integrals 
involving / M (m 1 ,m2;m 3 ), left to be carried out in Eqs. ( |C.9| ), (|C.10|) . These integrals 
are of course well-defined and finite. However, if m 3 < \mi — m 2 | or m 3 > m x + m 2 , 
they involve integration over logarithmic singularities. In our numerics, we found that 
the integration is more effective if we factorise out the singularities explicitly. Suppose 
we, for instance, first perform the integral over q. If mi 7^ 0, we write 

(m 3 -m x -m 2 ), 



2m? 



i-|(j» 2 + m\) 1,2 (ml - 2ral(m\ + m\) + (ml 
(q — a — b)(q + b — a) 



2m{ 

f p , q (m 1 ,m 2 ;m 3 ,) 



m 9 



,2\2 



1/2 



In 



(q + a + b)(q + a - b) 
and there are then singularities at q — \a + b\, \a — b\. If mi 
Am\p 2 — (m 2 — ml) 2 



Ap(m 2 



ml) 



f p , q (0,m 2 ;m 3 ) 



In 



0, we write 

q-io 



and there is a singularity at q = \qo\. If m 2 = m 3 , f p , q (0, m; m) 



q + Qo 
= 0. 



(C.12) 
(C.13) 
(C.14) 

(C.15) 



One heavy, one light mass. Let us now consider in more detail some special cases of 
Hb(mi, m 2 , m 3 ) needed in the main part of this paper. For the consideration in Sec. 0, 
we need to know how H^mn, mu, ttiq) behaves for small m#, m\j. We claim that 
there is a linear term oc mu (modulo logarithms). Since the result is symmetric 
in mn, "rn v , it is enough to consider mu <C mq. 

Non-analytic terms can only arise from a zero Matsubara mode. Thus, 

rf 3-2e„ 



Hi 



V 



1 



1 



(27r) 3 - 2e p 2 + m 2 H %ql + (p + q) 2 + mfj q 2 + q 2 + m 2 Q 



(C.16) 



Let us denote the latter integral by II(p). It is then obvious that the leading behaviour 
must be 



m H <^m Q 



d 



3-2e 



V 



(2ir) 3 ~ 2 * p 2 + m 2 H 



n(o) = I 3d {m H )- 



m f 



my 



mf 



hd(m H ) -I{m Q ) + hd(mu) + 4.^o(0) 



I(m v ) - I(m Q ) 



1 2 
'mfj mfj 



^m, 



2 ■ 

Q 



J>2 



(C.17) 
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Figure 3: The finite part of H b (mH,mQ,0) — H b (0,rriQ,0) (the zero temperature con- 
tribution has been subtracted) with solid lines, compared with the "linear term" in 
Eq. (|C.17| ) with dashed lines. The results have been divided by T 2 /(An) 2 . 



Indeed, the remainder, 

~/ (£^^ + U [ n(p) " n(0) ]- (ai8) 

behaves at small p as / dpp A / (p 2 + m 2 H ). This is IR finite even after expanding in m 2 H , 
therefore there cannot be any further linear contributions. 

In order to verify this behaviour explicitly, we set m\j = and compare the finite 
part of Eq. (|C.17|) with a numerical evaluation of the finite part of H b (mH,mQ,0) in 
Eq. (|C.9|) . We fix p, = T, and choose tuq/T = 1.0,2.0,4.0. The result is shown in 
Fig. Note that in Eq. ( |C.17|) , one must include a contribution ~ m H \nm H , arising 
when the (9(e) part of I 3d (m H ) (cf. Eq. ( |A.12| )) combines with the 1/e pole in /(m<g). 
From the perfect agreement at small m^/T in Fig. [|, we conclude that for m#/T 1 
the behaviour is indeed according to Eq. ( |U.17| ). 



Two equal heavy masses. Finally, let us consider the special case needed in Sec. ffl 
H b (m,m,0) (cf. Eq. ( |4.2|) ). The T = part in Eq. ( |C.11| ), related to 2-loop vacuum 
renormalization of m 2 H (fl), is 



iy vac (m,m,0) = - 




(C.19) 
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H 



finite 



1.0 - 




no expansion 



\ low-T expansion 
\ 

v«- high-T expansion 

\ i I 



0.0 



2.0 
m/T 



4.0 



1.0- 



0.0 



low-T expansion/no expansion 



high-T expansion/no expansion 



0.0 



2.0 
m/T 



4.0 



Figure 4: Left: The function Tt = H &nite in Eq. ( |C.21| ), along with the high-T and 
low-T limits in Eq. ( |C.22| ). Right: The ratio of the high-T and low-T limits to the 
non-expanded result. 



Using Eq. ( |A.6| ) as well as the simple expressions for n(0; m 2 , m 2 ), II(— m 2 ; m 2 , 0) ob- 
tained from Eqs. dCTTf) , flCH ), we then get from Eq. (|09l) 



H b (m, to, 0) = H vac (m, m, 0) + /i 



-4c 



T 2 



(4tt) 



1 _ /m 
+ 2i 



12 



•T' V T' 6 
where X\ , X 2 are from Eq. ( |A.8|) , and 

2 



T 

+(4 _ 21n2)Il Q_ l2 Q + i( 7E _ ln2 _^|) +<) 



(C.20) 



= 4 HpP [% q n T p}n T 9 h n P -±«. (C.21) 
v T ' n 2 Jo Jo oj p uj q p — q 

The function TC(m/T), numerically very easily evaluated, is plotted in Fig. f|, together 
with a comparison with the limiting values 



T "3 +7s 



C(2)7 : 



l 

— ( 

2vr 



-2m/T 



(C.22) 
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